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We describe a standard star catalog constructed using multiple SDSS photo- 
metric observations (at least four per band, with a median of ten) in the ugriz 
system. The catalog includes 1.01 million non- variable unresolved objects from 
the equatorial stripe 82 (|<5/2ooo| < 1-266°) in the RA range 20h 34m to 4h 00m, 
and with the corresponding r band (approximately Johnson V band) magnitudes 
in the range 14-22. The distributions of measurements for individual sources 
demonstrate that the photometric pipeline correctly estimates random photo- 
metric errors, which are below 0.01 mag for stars brighter than (19.5, 20.5, 20.5, 
20, 18.5) in ugriz, respectively (about twice as good as for individual SDSS 
runs). Several independent tests of the internal consistency suggest that the spa- 
tial variation of photometric zeropoints is not larger than ~0.01 mag (rms). In 
addition to being the largest available dataset with optical photometry internally 
consistent at the ~1% level, this catalog provides practical definition of the SDSS 
photometric system. Using this catalog, we show that photometric zeropoints for 
SDSS observing runs can be calibrated within nominal uncertainty of 2% even 
for data obtained through 1 mag thick clouds, and demonstrate the existence of 
He and H white dwarf sequences using photometric data alone. Based on the 
properties of this catalog, we conclude that upcoming large-scale optical surveys 
such as the Large Synoptic Survey Telescope will be capable of delivering robust 
1% photometry for billions of sources. 

Subject headings: surveys — catalogs — standards — methods: data analysis - 
techniques: photometric — instrumentation: photometers 

1. Introduction 

Astronomical optical photometric data are usually calibrated using sets of standard stars 
whose brightness is known from previous work. The most notable modern optical standard 
star catalogs are Landolt standards (Landolt 1992) and Stetson standards (Stetson 2000, 
2005). Both are reported on the Johnson-Kron-Cousins system (Landolt 1983 and references 
therein). The Landolt catalog provides magnitudes accurate to 1-2% in the UBVRI bands 
for ~500 stars in the V magnitude range 11.5-16. Stetson has extended Landolt's work to 
fainter magnitudes, and provided the community with ~l-2% accurate magnitudes in the 
BVRI bands for ~15,000 stars in the magnitude range V < 20. Most stars from both sets are 
distributed along the Celestial Equator, which facilitates their use from both hemispheres. 

The data obtained by the Sloan Digital Sky Survey (SDSS; York et al. 2000) can be 
used to extend the work by Landolt and Stetson to even fainter levels, and to increase 
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the number of standard stars to over a million. In addition, SDSS has designed its own 
photometric system (ugriz, Fukugita et al. 1996) which is now in use at many observatories 
worldwide. This widespread use of the ugriz photometric system motivates the construction 
of a large standard star catalog with ~1% accuracy. As a part of its imaging survey, SDSS 
has obtained many scans in the so-called Stripe 82 region, which is defined by | £72000 1 < 
1.266° and RA approximately in the range 20h - 4h. These repeated observations can be 
averaged to produce more accurate photometry than the nominal 2% single-scan accuracy 
(Ivezic et al. 2004a). 

The catalog and methods presented here have some similarity with an effort by Pad- 
manabhan et al. (2007), who developed a new calibration algorithm that simultaneously 
solves for the calibration parameters and relative stellar fluxes using overlapping SDSS ob- 
servations (the so-called iibercalibration method). The algorithm decouples the problem of 
"relative" calibrations (i.e. producing an internally consistent system), from that of "ab- 
solute" calibrations (i.e. typing the internal system to a physical flux scale): the absolute 
calibration is reduced to determining a few numbers for the entire survey. Here we also 
decouple "relative" and "absolute" calibrations and use overlapping observations. The main 
difference between their work and this paper is that they are concerned about calibrating 
the entire SDSS survey (~8500 deg 2 ; most of the surveyed area has at most two overlapping 
observations), while we concentrate here on a much smaller area (~300 deg 2 ) with an aver- 
age of ten overlapping observations. We also determine flatfield corrections (relative to the 
"standard" survey reductions) using different methods: Padmanabhan et al. minimize er- 
rors in relative photometry of multiply observed stars, while we require that the stellar locus 
remains fixed in multi-dimensional color space. An advantage of the catalog presented here 
is better averaging of various photometric errors thanks to a larger number of observations, 
which comes at the expense of a much smaller cataloged area. It is encouraging that the 
results of these two complementary approaches agree in the regions of sky common to both 
catalogs at the claimed level of accuracy (~1%). 

Additional motivation for the analysis of repeated scans and their impact on photometric 
accuracy comes from upcoming large-scale optical surveys such as the Dark Energy Survey 
(Flaugher et al. 2007), Pan-STARRS (Kaiser 2002) and the Large Synoptic Survey Telescope 
(Tyson 2002, LSST hereafter). For example, the LSST science requirements document^ calls 
for a photometric system that is internally consistent across the sky at the 1% level. The 
SDSS Stripe 82 repeated scans can be used to gauge the plausibility of delivering such a 
system. 



1 Available from |http://www.lsst.org/Science/lsst_baseline.shtml 
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We describe the construction and testing of a standard star catalog in §2, and illustrate 
its several use cases in 53. We discuss our results in 54. 



2. The Construction of SDSS Stripe 82 Standard Star Catalog 

2.1. Overview of SDSS imaging data 

SDSS is using a dedicated 2.5m telescope (Gunn et al. 2006) to provide homogeneous 
and deep (r < 22.5) photometry in five bandpasses (Fukugita et al. 1996; Gunn et al. 1998, 
2006; Smith et al. 2002; Hogg et al. 2002) repeatable to 0.02 mag (root-mean-square scatter, 
hereafter rms, for sources not limited by photon statistics, Ivezic et al. 2003) and with a 
zeropoint uncertainty of ~0. 02-0. 03 (Ivezic et al. 2004a). The survey sky coverage of close 
to ~10,000 deg 2 in the Northern Galactic Cap, and ~300 deg 2 in the Southern Galactic 
Hemisphere, will result in photometric measurements for well over 100 million stars and a 
similar number of galaxieaj. Astrometric positions are accurate to better than 0.1 arcsec 
per coordinate (rms) for sources with r < 20. 5 m (Pier et al. 2003), and the morphological 
information from the images allows reliable star-galaxy separation to r ~ 21.5 m (Lupton et 
al. 2002, Scranton et al. 2002). 

Data from the imaging camera (thirty photometric, twelve astrometric, and two focus 
CCDs, Gunn et al. 1998) are collected in drift scan mode. The images that correspond to the 
same sky location in each of the five photometric bandpasses (these five images are collected 
over ~5 minutes, with 54 sec for each exposure) are grouped together for simultaneous 
processing as a field. A field is defined as a 36 seconds (1361 pixels, or 9 arcmin, see 
Stoughton et al. 2002) long and 2048 pixels wide (13 arcmin) stretch of drift-scanning data 
from a single column of CCDs (sometimes called a scanline, for more details please see 
Stoughton et al. 2002; Abazajian et al. 2003, 2004, 2005; Adelman-McCarthy et al. 2006). 
Each of the six scanlines (called together a strip) is 13 arcmin wide. The twelve interleaved 
scanlines (or two strips) are called a stripe (~2.5° wide). 

2.2. The photometric calibration of SDSS imaging data 

SDSS 2.5m imaging data are photometrically calibrated using a network of calibration 
stars obtained in 1520 41.5x41.5 arcmin 2 transfer fields, called secondary patches. These 



2 The recent Data Release 5 (Adelman-McCarthy et al. 2007, in prep) lists photometric data for 215 
million unique objects observed in 8000 deg 2 of sky; see http://www.sdss.org/dr5/! 
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patches are positioned throughout the survey area and are calibrated using a primary stan- 
dard star network of 158 stars distributed around the Northern sky (so-called USNO stan- 
dards, Smith et al. 2002). The primary standard star network is tied to an absolute flux 
system by the single F0 subdwarf star BD+17°4708, whose absolute fluxes in SDSS filters 
are taken from Fukugita et al. (1996). The secondary patches are grouped into sets of four, 
and are observed by the Photometric Telescope (hereafter PT; Tucker et al. 2006) in parallel 
with observations of the primary standards. A set of four patches spans all 12 scanlines of 
a survey stripe along the width of the stripe, and the sets are spaced along the length of 
a stripe at roughly 15 degree intervals, which corresponds to an hour of scanning at the 
sidereal rate. 

SDSS 2.5m magnitudes are reported on the "natural system" of the 2.5m telescope 
defined by the photon-weighted effective wavelengths of each combination of SDSS filter, 
CCD response, telescope transmission, and atmospheric transmission at a reference airmass 
of 1.3 as measured at APCfl The magnitudes are referred to as the ugriz system (which 
differs from the "primed" system, u'g'r'i'z', that is defined by the USNO standards^). The 
reported magnitude^ are corrected for the atmospheric extinction (using simultaneous ob- 
servations of standard stars by the PT) and thus correspond to measurements at the top of 
the atmosphere^] (except for the fact that the atmosphere has an impact on the wavelength 
dependence of the photometric system response). The magnitudes are reported on the AB 
system (Oke & Gunn 1983) defined such that an object with a specific flux of F^=3631 Jy 
has m = (i.e. an object with F i/ =const. has an AB magnitude equal to the Johnson V 
magnitude at all wavelengths). In summary, given a specific flux of an object at the top of 
the atmosphere, F V (X), the reported SDSS 2.5m magnitude in a given band, b=(u, g,r,i, z), 
corresponds to (modulo random and systematic errors, which will be discussed later) 



Transmission curves for the SDSS 2.5m photometric system are available at 
http: / / www.sdss.org/ dr5 / instruments / imager 



4 For subtle effects that led to this distinction, please see Stoughton et al. (2002), Smith et al. (2002), 
and http: / /www. sdss.org/ dr5 / algorithms/nuxcal.html 

5 SDSS uses a modified magnitude system (Lupton, Szalay & Gunn 1999), which is virtually identical to 
the standard astronomical Pogson magnitude system at high signal-to-noise ratios relevant here. 

6 The same atmospheric extinction correction is applied irrespective of the source color; the systematic 
errors this introduces are probably less than 1% for all but objects of the most extreme colors. 



- 6- 



where 

r F„(A)&(A)d\. (2) 



Here, </>&(A) is the normalized system response for the given band, 

with the overall atmosphere + system throughput, S&(A), available from the website given 
above (0&(A) for the SDSS system are shown in Figure El see also § 12.5.21) . We reiterate 
that the normalization of reported magnitudes corresponds to a source at the top of the 
atmosphere, while the throughput </>&(A) includes the transmission of a standard atmosphere 
at a fiducial airmass of 1.3. Note also that it is only the shape of S&(A), and not its overall 
normalization, that needs to be known to compute expected SDSS magnitudes of a source 
with given F U (X). That is, the SDSS photometric system is fully defined by the five dimen- 
sionless functions 0&(A) (by definition, J (p b dX = 1, see eq. [3]). In reality, for each ugriz band 
there are six devices in the SDSS camera (Gunn et al. 1998) whose 0& are slightly different 
(see §1232]). 

The quality of SDSS photometry stands out among available large-area optical sky 
surveys (Ivezic et al. 2003, 2004a; Sesar et al. 2006). Nevertheless, the achieved accuracy 
is occasionally worse than the nominal 0.02-0.03 mag (root-mean-square scatter for sources 
not limited by photon statistics). Typical causes of substandard photometry include an 
incorrectly modeled point spread function (PSF; usually due to fast variations of atmospheric 
seeing, or lack of a sufficient number of the isolated bright stars needed for modeling the 
PSF), unrecognized changes in atmospheric transparency, errors in photometric zeropoint 
calibration, effects of crowded fields at low Galactic latitudes, an undersampled PSF in 
excellent seeing conditions (< 0.8 arcsec; the pixel size is 0.4 arcsec), incorrect flatfield or 
bias vectors, scattered light erroneously included in flatfield, etc. Such effects can conspire to 
increase the photometric errors to levels as high as 0.05 mag (with a frequency, at that error 
level, of roughly one field per thousand). However, when multiple scans of the same sky region 
are available, many of these errors can be minimized by properly averaging photometric 
measurements. 



2.3. The Choice of Cataloged Magnitudes 



The SDSS photometric pipeline (photo, Lupton et al. 2002) measures several types 
of magnitudes, including aperture, PSF, and model magnitudes. Here we briefly describe 
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each type of magnitude (for more details see Stoughton et al. 2002, and the SDSS website 
www.sdss.org) and justify the choice of PSF magnitudes for catalog construction. 

2.3.1. Aperture magnitudes 

Aperture magnitudes computed by photo are based on the flux contained within the 
aperture with a radius of 7.43 arcsec. While an aperture magnitude is the most robust flux 
estimate at the bright end (because it is essentially seeing independent), these magnitudes do 
not have good noise properties at the faint end where sky noise dominates (e.g. for a given 
maximum photometric error, PSF magnitudes reach 1-1.5 mag fainter than do aperture 
magnitudes). In order to improve the depth of the standard star catalog, we opt not to use 
aperture magnitudes, except for quality tests at the bright end. 

2.3.2. Point spread function magnitudes 

The point spread function (PSF) flux is computed using the PSF as a weighting function. 
While this flux is optimal for faint point sources (in particular, it is vastly superior to aperture 
photometry at the faint end), it is also sensitive to inaccurate PSF modeling as a function of 
position and time. Even in the absence of atmospheric variations, the SDSS telescope and 
camera optics deliver images whose FWHMs vary by up to 15% from one side of a CCD to 
the other; the worst effects are seen in the chips farthest from the optical axis. Moreover, 
since the atmospheric seeing varies with time, the delivered image quality is a complex two- 
dimensional function even on the scale of a single frame. Without accounting for this spatial 
variation, the PSF photometry would have errors up to 0.10-0.15 mag. The description of 
the PSF is also critical for star-galaxy separation and for unbiased measures of the shapes 
of nonstellar objects. 

The SDSS imaging PSF is modeled heuristically in each band and each camera column 
using a Karhunen-Loeve (KL) transform (Lupton et al. 2002). Using stars brighter than 
roughly 20 th magnitude, the PSF from a series of five frames is expanded into eigenimages and 
the first three terms are retained. The variation of these coefficients is then fit up to a second 
order polynomial in each chip coordinate. The failure of this KL expansion, typically due to 
insufficient number of stars, or exceedingly complex PSF, results in occasional problems with 
PSF photometry. The main failure mode is inaccurate determination of aperture corrections 
which statistically tie PSF magnitudes to aperture magnitudes using bright stars. 
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2.3.3. Model magnitudes 

Just as the PSF magnitudes are optimal measures of the fluxes of stars, the optimal 
measure of the flux of a galaxy uses a matched galaxy model. With this in mind, the 
photometric pipeline fits two models to the two-dimensional image of each object in each 
band: a pure deVaucouleurs profile and a pure exponential profil^E Because the models 
are convolved with a double-Gaussian fit to the PSF, the seeing effects are accounted for. 
Aperture corrections are applied to make these model magnitudes equal the PSF magnitudes 
in the case of an unresolved object. 



2.3.4- The choice of magnitudes for the standard star catalog 

A comparison between aperture, PSF and model magnitudes for unresolved sources is 
done automatically for every SDSS observing run (runQA pipeline, Ivezic et al. 2004a). 
An analysis of over 200 runs indicate that model magnitudes are more robust than PSF 
magnitudes: PSF magnitudes show systematic offsets of 0.05 mag from aperture magnitudes 
three times more often than do model magnitudes (roughly once per thousand fields), as 
model fits have more degrees of freedom. On the other hand, an analysis of repeated scans 
indicates that estimates of photometric errors by the photometric pipeline are more accurate 
for PSF magnitudes (agreeing at the 10% level with the measured values, see Ivezic et al. 
2003; Scranton et al. 2005) than for model magnitudes (which are smaller than the measured 
values by typically 30-50%). Because the rejection of likely variable sources, which relies on 
accurate photometric error estimates, is an important step in the construction of the standard 
star catalog (see below), we choose to use PSF magnitudes to construct the catalog. 



2.4. Catalog Construction 

Using 58 SDSS-I runs from stripe 82 (approximately 20h < ctj 2 ooo < 04h and |5j2ooo| < 
1.266°, but not all runs extend over the entire right ascension range) obtained in mostly 
photometric conditions (as indicated by the PT calibration residuals, infrared cloud camera^, 
and tests performed by the runQA pipeline), candidate standard stars from each run are 
selected by requiring 



7 For more details see |http: / / www.sdss.org/dr5 /algorithms/photometry.html 



'For more details about the camera see http://hoggpt.apo.nmsu.edu/irsc/irsc_doc and Hogg et al. (2002). 
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1. that objects are classified as STAR (based on the difference between model and PSF 
magnitudes); this morphological classification really means unresolved (point) sources 
(e.g. quasars are also included), 

2. that they have quoted photometric errors in the PSF magnitude (as computed by the 
photometric pipeline) smaller than 0.05 mag in at least one band, and 

3. that the processing flags BRIGHT, SATUR, BLENDED, or EDGE are not set in any 
bancB 

These criteria select unsaturated point sources with sufficiently high signal-to-noise ratio per 
single observation to approach final photometric errors of 0.02 mag or smaller. 

After matching all detections of a single source across runs (using a 1 arcsec matching 
radius), various photometric statistics such as unweighted mean, median, their standard 
errors, root-mean-square scatter, number of observations, and \ 2 P er degree of freedom, are 
computed for magnitudes in each band. We use errors reported by the photometric pipeline 
to compute \ 2 an d note that systematic errors common to all runs do not contribute to its 
value. This initial catalog of multi-epoch observations includes 1.4 million point sources with 
at least 4 observations in each of the g, r and i bands. The median number of observations 
per source and band is 10, and the total number of photometric measurements is ~57 million. 

The distributions of the median magnitudes, their standard errors, x 2 an d the number 
of observations for a subset of these sources are shown in Figure [TJ The random errors 
in the median magnitude (computed as 0.928*IQR/\M^ — 1, where IQR is the 25%- 75% 
interquartile range of the individual measurement distribution and N is the number of mea- 
surements; note that the error of the median for a Gaussian distribution is 25% larger than 
the error of the mean, Lupton 1993) are below 0.01 mag at the bright end. These errors are 
reliably computed by the photometric pipeline, as indicated by the \ 2 distributions. The 
distributions of these sources in color-magnitude and color-color diagrams, constructed using 
median magnitudes, are shown in Figure [2j For a detailed interpretation of these diagrams, 
please see Lenz et al. (1998), Fan (1999), Finlator et al. (2000), Helmi et al. (2003), and 
Ivezic et al. 2006. As evident, the sample is dominated by stars. 



9 For more details about processing flags see 
Stoughton et al. (2002). 



http: / /www.sdss.org/ dr5 /products / catalogs / flags.html 



and 
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2.4-1- Selection of Candidate Standard Stars 

Adopted candidate standard stars must have at least 4 observations and, to avoid vari- 
able sources, x 2 l ess than 3 in the gri bands (the same requirements are later applied in 
the u and z bands when using the catalog for calibration, as we discuss further below). We 
also limit the right ascension to the range from 20h 34m to 4h 00m, which provides a simple 
areal definition (together with |5j2ooo| < 1-266°) of a 282 deg 2 large rectangular region, while 
excluding only a negligible fraction of stars. With the final condition that the standard error 
for the mean magnitude in the r band is smaller than 0.05 mag, these requirements result 
in a catalog with slightly over a million sources (1,006,849) with 42 million photometric 
measurements. Of those, 563,908 have a random error for the median magnitude in the r 
band smaller than 0.01 mag, and for 405,954 stars this is true in all three of the gri bands. 
Subsets of 92,905 and 290,299 stars satisfy these requirements in the ugri and griz bands, 
and 91,853 stars satisfy this in all five bands. The distributions of candidate standard stars 
that satisfy the above selection criteria in all five bands in color-magnitude and color-color 
diagrams are shown in Figure 

For comparison, the distribution of sources that were rejected as variable (x 2 greater 
than 3 in at least one of the gri bands) in color-magnitude and color-color diagrams is 
shown in Figure HI As evident from a comparison with Figure [31 the distribution of variable 
sources in color-color diagrams is markedly different from that of non-variable sources. It is 
especially striking how low-redshift (z < 2.2) quasars are easily detected by their variability 
(for more details see Ivezic et al. 2004b). However, it is fairly certain that not all variable 
sources are recognized as such because of the limited number of repeated observations (~10). 
For example, an eclipsing binary with much shorter eclipse duration than the orbital period 
could easily escape detection. Analysis of the variable subsample is presented in a companion 
paper (Sesar et al. 2007, in prep.). 

The sky density of all the sources and those selected as non-variable are shown in 
Figure El At high Galactic latitudes ~ 60°) the fraction of candidate standard stars 
sources is ~80%. 

2.5. Systematic Photometric Errors 

Photometric errors computed by the photometric pipeline provide a good estimate of 
random errors in SDSS photometry, as demonstrated by the x 2 distributions shown in Fig- 
ure [H However, the measurements are also subject to systematic errors such as spatial 
dependence of the internal zeropoints (calibration errors), and the overall deviations of the 
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internal SDSS zeropoints from an AB magnitude scale. Formally, the true AB magnitude of 
an object (defined by eq{T]) in a given band, m true , can be expressed as 

m t rue = m ca t + S m (RA, Dec) + A m , (4) 

where m cat is the cataloged magnitude, <5 m (RA,Dec) describes the spatial variation of the 
internal zeropoint error around A m (thus the average of 5 m over the cataloged area is by 
construction), and A m is the overall (spatially independent) deviation of the internal SDSS 
system from a perfect AB system (the five values A m are equal for all the cataloged objects). 
Here we ignore systematic effects, e.g., device non-linearity and bandpass variations between 
different camera columns, which depend on individual source properties such as brightness 
and colors (but see §2.5.21 below). 

The spatial variation of the internal zeropoint error can be separated into color errors, 
relative to a fiducial band, say r, and an overall "gray" error (e.g. unrecognized temporal 
changes in atmospheric transparency due to gray clouds) 

<5 m (RA, Dec) = 5 r (RA, Dec) + S mr (RA, Dec). (5) 

Below we discuss methods for estimating both the "gray" error <5 r (RA,Dec) and the color 
errors <5 mr (RA, Dec). 

The deviation of the internal SDSS system from a perfect AB system, A m , can also be 
expressed relative to the fiducial r band 

A m A r -|- A mr . (6) 

The motivation for this separation is twofold. First, A mr can be constrained by considering 
the colors (spectral energy distributions) of objects, independently from the overall flux scale 
(this can be done using both external observations and models). Second, it is difficult to find 
a science result that crucially depends on knowing the "gray scale" offset, A r , at the 1-2% 
level. On the other hand, knowing the "band-to-band" offsets, A mr , with such an accuracy is 
important for many applications (e.g., photometric redshifts of galaxies, type la supernovae 
cosmology, testing of stellar and galaxy models). 

Fitting SDSS spectra of hot white dwarfs to models, Eisenstein et al. (2006) determined 
the AB color corrections A mr to be —0.04, 0.00, —0.015 and —0.015 mag for m = ugiz, 
respectively, with an uncertainty of ~0. 01-0. 02 mag. It may be possible to determine these 
corrections with an uncertainty of ~0.01 mag, and such efforts are in progress (J. Marriner, 
priv. comm.). 

The overall "gray" flux scale calibration error, A r , is determined by the accuracy of 
the absolute flux calibration of fundamental standard BD+17°4708 (Fukugita et al. 1996), 
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the accuracy of tying the primary standard star network to BD+17°4708, the accuracy of 
transfering the primary standard star network to the secondary standard star network, and 
the accuracy of the calibration of the survey imaging data using the secondary standard star 
network. Given these numerous sources of error, it seems unlikely that A r < 0.02 mag. On 
the other hand, formal analysis of all the error contributions, as well as a direct comparison 
to HST observations of hot white dwarfs (Eisenstein et al. 2006), suggest that A r does not 
exceed 0.05 mag. Note, however, that all these uncertainties in the definition and transfer 
of the standard star network become moot if one accepts that 

1. A r does not need to be known exquisitely well for most scientific applications. Even if 
it does, this is just a single number that modifies the cataloged photometry for all the 
sources and all the bands in the same fashion. 

2. Uncertainties in the determination of A mr are of the order 0.01 mag. 

3. <5 m (RA, Dec) can be constrained, or corrected for, at the 0.01 mag level. 

In other words, the band-to-band calibration can be fixed by adopting A mr , determining and 
correcting for <5 m (RA, Dec) guarantees internal consistency, and the only remaining relatively 
free parameter is A r . Such a system is then no longer defined by a set of celestial standards 
but rather by the functions <pb- The catalog presented here is one realization of such a 
photometric system. 



2.5.1. Determination of 5 m (RA, Dec) 

We now proceed to describe methods for constraining 5 m (RA, Dec). The region covered 
by the SDSS Stripe 82 is an elongated rectangle with an aspect ratio of 1:50, and with the 
long side parallel to the Celestial Equator. Because of this large aspect ratio, and because 
different effects contribute to the RA and Dec dependences of 5 m , we assume that it can be 
expressed as a sum of independent functions of RA or Dec, respectively 

URA, Dec) = <S#(Dec) + <^(RA), (7) 

with (<^(Dec)) RA = and (<^f (RA)) Dec = 0, where (6) x denotes the average of 5 over 
direction x. 

The first term, 8^ (Dec), is dominated by the errors in the flatfield vectors (for drift 
scanning, flatfield corrections are one-dimensional). The flatfield determination for SDSS 
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was difficult due to scattered light! 10 l and there are probably systematic errors in the stellar 
photometry at the 0.01 mag level in the griz bands and the 0.02 mag level in the u band 
(perhaps somewhat larger at the edge of the imaging camera), as we will demonstrate below. 
Since these systematic errors do not cancel when averaging many observing runs (because 
most stars are always observed by the same CCD and fall on roughly the same position 
within the CCD), 8^ (Dec) could be as large as ~ 0.01 — 0.02 mag on spatial scales much 
smaller than the chip width of 13 arcmin. 

The second term, 5^{RA), is dominated by unrecognized fast variations of atmospheric 
extinction (e.g. due to cirrus), because for each observing run only a single zeropoint per 
CCD is determinecQ While such variations are uncorrelated for different runs, it is possible 
that they do not average out fully at the 1% level. 



2.5.2. The six ugriz photometric systems 

Before we describe methods for determining 5 m (RA, Dec), we address the bandpass 
differences between the six camera columns. The bandpasses had been measured for each 
CCD using a monochromator (M. Doi et al. 2007, in prep.) and found not to be identical, 
as shown in Figure El These differences between the bandpasses induce color term errors in 
the reported SDSS photometry because the magnitudes of calibration stars obtained by the 
PT are transformed to the 2.5m system using a single set of color termini In other words, 
the color difference between a blue star and a red star depends on which camera column the 
stars fall on: at the ~1% accuracy level discussed here, there are six SDSS ugriz systems. 

We used the measured response curves to generate synthetic ugriz photometry corre- 
sponding to six camera columns (via eq{T]) for 175 stars from the Gunn-Stryker atlas (1983). 
The differences between the predicted magnitudes for each camera column and the values 
generated with the response curves which define the SDSS systenj^l represent photometric 
corrections due to bandpass differences. These corrections may be as large as 0.02 mag for 



10 For details please see http://www.sdss.org/dr5/algoritlims/llatfield.litnil and Stoughton et al. (2002). 

11 Although not relevant for equatorial runs discussed here, slow change of atmospheric extinction due to 
varying airmass is accounted for; the response of the telescope and CCDs is stable at the <1% level on 
single-night timescales. 

12 See http: / /www. sdss.org/dr5 / algorithms/jeg_photometric_eq_drl .html 



13 Transmission curves for the SDSS 2.5m photometric system are available 



at 



http://www.sdss.org/dr5/instruments/imager 



-14- 



the reddest stars, but fortunately admit simple linear fits as a function of color (u — g for the 
u band, g — r for g and r, and r — i for i and z). The u band is the only one where piecewise 
fits are required (in the range 0.7 < u — g <1.3, the g — r color provides a better fit than the 
u — g color). We have applied these best-fit corrections to all sources in the catalog. The 
median and rms scatter for the distributions of corrections evaluated for all stars, and for 
each color and camera column combination, are listed in Table 1. 

The bandpass differences have the largest impact on the i — z color of red sources. The 
color term errors result in a rotation of the stellar locus in the i — z vs. r — i color-color 
diagram (see the bottom right panel in Figure [2]), and we utilize this fact to demonstrate 
the improvement in photometry due to applied corrections. We use the mean position of 
the stellar locus to "predict" the i — z color from the measured r — i color, and compute the 
difference between predicted and measured i — z colors separately for blue (0.1< r — i < 0.2) 
and red (0.8< r — i <1.4) stars, and in small bins of declination (cross-scan direction). 
The difference of these residuals (blue vs. red stars) effectively measures the locus position 
angle and is not sensitive to photometric zeropoint errors and flatfield errors (which can only 
induce locus shifts that have no effect on this test because they cancel out). 

The top panel in Figure [7] shows the median i — z residuals before applying corrections 
for different transmission curves, and the bottom panel shows results based on corrected 
photometry. The rms scatter decreases from 9 millimag to 3 millimag after applying the 
corrections. The remaining deviations of residuals from zero could be due to the fact that 
we have not corrected for dependence of the flatfields on source color (measurements of this 
dependence are not available). Because this is currently the only available SDSS catalog with 
photometry corrected for bandpass differences, it effectively represents a practical definition 
of the SDSS photometric system (i.e. it provides photometry for real sources on the sky; the 
system is formally defined by of course). 



2. 5. 3. Other Sources of Systematic Errors 

The non-linearity of the detectors (as a function of source brightness) has also been 
measured in situ by Doi et al. and found to be at most a couple of percent effect over the 
relevant dynamic rang These corrections are determined with a sufficient accuracy (<5 
millimag impact on photometry) and are already implemented in the photometric pipeline 
(i.e. before performing photometric calibration). 



For details see |http: / / www.sdss.org/dr5 /instruments/imager /nonlinearity.html 
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Similarly to the camera column-to-column bandpass differences, variations in the wave- 
length dependence of atmospheric transmissivity can also induce systematic errors that de- 
pend on source colors. Assuming a standard atmosphere, we find using synthetic photometry 
for stars from the Gunn-Stryker atlas that this effect can induce offsets of up to ~0.01 mag 
for the u — g and g — r colors when airmass is varied by 0.3 from its fiducial value of 1.3. How- 
ever, because all stripe 82 data are obtained at the same airmass, this effect is not relevant 
for the catalog discussed here. At least in principle, similar, and potentially larger, errors 
could be induced even at a constant airmass if the wavelength dependence of atmospheric 
transmissivity is significantly different from the assumed standard atmosphere. Given that 
such errors would probably average out, and that there are no available measurements of 
the wavelength dependence of atmospheric transmissivitjQ for SDSS data considered here, 
we ignore this effect hereafter. 



2.6. Determination of Flatfield Corrections 

We use two methods based on SDSS data to constrain <5 m (RA, Dec): a direct comparison 
with the secondary standard star network ( §2.6.2[) and a method based on stellar colors 
( §2.6.ip . Each of these methods has its advantages and disadvantages. We perform the tests 
of the final catalog quality using a method based on the photometric redshift relation for 
galaxies ( §2.7.2p . and also compare the SDSS photometry to an independent set of standards 
provided by Stetson (2000, 2005) in §EZ! 

A determination of the 5 m (RA, Dec) error from a direct comparison with the secondary 
standard star network (hereafter PT comparison) might be considered as the best method 
a priori However, it is quite possible that the secondary standard star network itself, off of 
which this catalog is calibrated, may induce a spatial variation of the photometric zeropoints 
at the 0.01 mag level (Smith et al. 2002). In addition, there are not enough stars to constrain 
8*1 (Dec) with sufficient spatial resolution (say, at least ~100 pixels, or ~0.01°). For example, 
there are ~20,000 secondary standards from Stripe 82 in the averaged catalog that are not 
saturated in the gri bands in the 2.5m scans, and have PT errors smaller than 0.03 mag 
(~8,000 stars are useable in the z band, and only ~3,000 in the u band). If these stars are 



15 Some handle on the stability of the wavelength dependence of atmospheric transmissivity can be obtained 
by studying first-order extinction coefficients determined by the photometric calibration pipeline. They show 
a cyclic variation during the year, with an rms scatter of residuals around the mean relation of 0.01 mag. 
The amplitude of the yearly variation of the r band extinction coefficient is ^20% about the mean value, and 
the wavelength dependence of the variation appears consistent with the addition of a gray opacity source 
during summer months. 
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binned in the declination direction every 0.01° (250 90-pixels wide bins), 5 f J in each bin can 
be constrained to about ~0.005 mag (0.01 mag in the u band). This is barely sufficient in the 
gri bands, and cannot provide satisfactory constraints on the flatfielding errors in the u band 
(where unfortunately these errors are the largest). Similarly, 5^*(RA) can be constrained 
in 0.5° wide right ascension bins with a similar accuracy, but the secondary standard stars 
are not uniformly distributed in right ascension. For these reasons, we combine the PT 
comparison with the stellar locus method to determine flatfield corrections. 



2.6.1. Color corrections from the stellar locus method 

The stellar distribution in color-color space at high Galactic latitudes^ (|fe| > 30) is 
remarkably uniform at the faint flux levels probed by SDSS, as discussed in detail by Ivezic 
et al. (2004a). Systematic photometric errors, other than an overall gray error, manifest 
themselves as shifts in the position of the stellar locus that can be tracked using the four 
principal colors (swxy) defined by Ivezic et al. (2004a). These colors are linear combinations 
of magnitudes, 

P 2 = Au + Bg + Cr + Di + Ez + F, (8) 

where P 2 = s, w, x, y, and measure the distance from the center of the locus in various two- 
dimensional projections of the four- dimensional stellar color distributions (s: perpendicular 
to the blue part of the locus in the g-r vs. u-g plane, w. perpendicular to the blue part in the 
r-i vs. g-r plane, x: perpendicular to the red part, with g — r ~ 1.4, in the r-i vs. g-r plane, 
and y: perpendicular to the locus in the i-z vs. r-i plane). The matrix of coefficients A — F 
is listed in Table 2 (for more details see Ivezic et al. 2004a). Of course, the measurements 
must be corrected for the effects of interstellar dust extinction; we use maps provided by 
Schlegel, Finkbeiner & Davis (1998, hereafter SFD98). The properties of two of these colors 
(s and w) are illustrated in Figure 

The fact that the median principal colors are close to zero shows that the averaging 
procedure did not induce any shifts in zeropoints compared to the average of 291 SDSS runs 
which were used to define the principal colors (Ivezic et al. 2004a). The same conclusion is 
reached by comparing the averaged photometry with the secondar y s tandard star network: 
the median photometric residuals at the so-called crossing colord_| are 4, 6, 3, 2 and 2 



16 At low Galactic latitudes several effects, discussed below, prevent the use of this method for calibration 
purposes. 

17 Crossing colors are roughly the median colors of the observed stellar population, for details see 
http : / / www. sdss . org/dr5/ algorithms /j eg_photometric_eq_dr 1 . html 
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millimag in the ugriz bands, respectively. Yet another test is a direct comparison of averaged 
photometry with single-epoch photometry. Using the SDSS Data Release 5 photometry, we 
find that the largest median magnitude difference between the two sets is 2 millimag in the 
u band. 

It is noteworthy that the widths of the principal color distributions (i.e. the thickness 
of the stellar locus) constructed with averaged photometry are much smaller than when 
using single-epoch data (see the bottom four panels in Figure ED- Indeed, all four principal 
color distributions are "resolved" using this high quality photometry (see the fifth column 
in Table 3). 

Because the intrinsic widths of the principal color distributions are so small, principal 
colors can be used to efficiently track local calibration problems using a small number of 
stars, allowing a high spatial resolution. That is, we require that the locus not move in the 
multi-dimensional color space. In practice, the deviations of the principal colors from can 
be inverted, using an appropriate closure relation (see the next section), to yield flatfield 
corrections (Ivezic et al. 2004a). With bins 0.01° wide in the declination direction, or 1° 
wide in the right ascension direction, the flatfield corrections can be determined with an 
accuracy of 5 millimag, or better. 

2.6.2. Gray corrections from the comparison with secondary standard star network 

The main advantage of the stellar locus method is that it can constrain 8 mr with high 
spatial resolution. However, it is insensitive to gray errors, parametrized by 5 r (e.g. an 
overall gradient of photometric zeropoints in the declination direction that is the same in all 
five bands would have no effect on stellar colors). On the other hand, the PT comparison can 
constrain 5 r , but it does not provide enough spatial resolution to derive flatfield corrections, 
especially in the u and z bands. Therefore, we combine these two methods to derive flatfield 
corrections 5^ (Dec). 

The median differences between the averaged 2.5m photometry and PT photometry for 
secondary standard stars in the gri bands are shown in the top panel in Figure [9j The 
median differences are computed for 0.01° wide bins, and then smoothed by a triangular 
filter (jji is replaced by 0.25 * {y%-\ + tyi + y%+i))- The residuals in all three bands display 
similar behavior and imply about 0.02 peak-to-peak variation between the center and edges 
on each CCD (resulting in about 6 millimag rms contribution to the overall errors), as well 
as an overall 0.01-0.02 mag tilt. These systematic errors may be due to imperfect flatfield 
vectors used to reduce the data, incorrectly determined scattered light correction (the two 
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are somewhat coupled in the data reduction procedures), or problems in PT itself (such as 
the PT flatfield). 

At face value, these residuals could be used to correct the averaged 2.5m photometry in 
each band (gri) separately. However, doing so introduces noise in stellar principal colors of 
about 5 millimag (rms) and suggests that the differences in photometric residuals between 
the three bands are dominated by PT measurement noise. On the other hand, the 2.5m vs. 
PT residuals do contain information about "gray" errors that cannot be determined using 
stellar locus. Hence, we take the mean value of the 2.5m vs. PT residuals in the gri bands 
to represent the 5 r flatfield correction, and apply it to the averaged 2.5m photometry in the 
r band. The applied r band correction is shown in the second panel in Figure [9] and has an 
rms scatter of 7 millimag (for 250 bins), with the largest correction less than 0.02 mag. 

In the second step, we use the stellar locus to derive the 5 mr corrections in each band 
(ugiz). The derivation of these corrections is essentially identical to the procedure described 
by Ivezic et al. (2004a). Also, together with a PT-based 5 r correction, this is essentially 
the same method as used to derive flatfield corrections for the whole SDSS surveyq. In 
particular, we used here the same closure relation (stellar locus method gives four equations 
for five unknowns), that is based on averaged 2.5m vs. PT residuals in the gri bands. The 
resulting flatfield corrections, 5 mr , in the ugiz bands are shown in the third and fourth panels 
in Figure [9] and summarized in Table 4. 

Due to low stellar counts and the strongest scattered light, the u band correction is 
expected to have the largest noise (~5-10 millimag), which is consistent with the observed 
behavior. It is thus likely that some of the variation on scales of ~0.01° is not real. On the 
other hand, it could be argued that systematic errors could actually be much larger on even 
smaller spatial scales, but get averaged out in 90 pixel wide bins. However, in addition to 
not having a reason to believe in such high spatial frequency effects (e.g. the sky background 
does not show any evidence for them), no additional scatter, except the expected statistical 
noise, is observed when the bin size is decreased by a factor of 4. 

Last but not least, it is important to emphasize that these corrections are not setting 
photometric zeropoints, but only correcting for variations in response across each CCD. As 
discussed above, the AB photometric zeropoints, relative to the fiducial r band, are effectively 
set by adopting values for A mr . 



For details see |http: / / www.sdss.org/dr5 /algorithms/ fiatfield.html 
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2.7. Tests of Catalog Quality 

2.7.1. The internal tests 

At least in principle, the same methods used to derive 6™ (Dec) could be used to derive 
<5^'(RA). However, in practice this is not possible for at least two reasons: first, the right 
ascension distribution of secondary standard stars is not as uniform as their declination 
distribution, and second, the assumption of the constancy of the stellar locus in color space 
is invalid along the "long" scan direction (discussed below). For these reasons, we only 
use the PT comparison and stellar locus methods to estimate the level of internal zeropoint 
variations with RA, and do not correct the data. In the next section, we will also use another 
method, based on galaxy colors, as an independent test of catalog integrity. 

Figure [TU] shows the median principal colors in bins of right ascension. As evident, the 
principal colors are close to zero for the right ascension in the range —25° < RA < 40°, 
but outside this range deviate significantly from zero. This does not necessarily indicate 
problems with photometric calibration, because the stellar locus method is expected to fail 
at low Galactic latitudes due to several reasons. First, the mean metallicity of stars increases 
at low Galactic latitudes and this change may affect the s and w colors. Second, at low 
latitudes red dwarfs are not behind the entire dust screen measured by the SFD98 maps (see 
J uric et al. 2006 for a discussion of this point), and thus the x color will be biased blue 
(i.e. the colors of red dwarfs are over-corrected for the ISM reddening). And third, at low 
latitudes the dust column increases fast (see Figure [TTT) and even small errors in the assumed 
wavelength dependence of the dust extinction, or the extinction itself as given by the SFD98 
maps, will have noticeable effects on principal colors. For these reasons, it seems plausible 
that the deviations seen in Figure [10] are not dominated by zeropoint errors. 

This conclusion is supported by the direct comparison of the averaged and PT photome- 
try (Figure [TT]) . For example, the largest median photometric residual between the averaged 
catalog and PT observations in the u band is ~0.02 mag (see Table 5), which is much smaller 
than the 0.1 mag discrepancy implied by the stellar locus method. 

Table 5 shows that the rms scatter of median photometric residuals (evaluated in 2° wide 
bins in the right ascension direction) between the averaged catalog and PT observations is 
<0.01 mag in all five bands. Some of that scatter must come from the PT data itself, and thus 
the true scatter of photometric zeropoints in the averaged catalog is probably even smaller 
than that listed in Table 5. In addition, Table 5 shows that the averaged catalog and PT 
measurements are on the same system to within a few millimags (using the recommended 
photometric transformations between the two telescopes listed at the SDSS website, see 
Section 
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We have also compared the catalog presented here to photometric reductions described 
by Padmanabhan et al. (2007). As discussed in the Introduction, they determined flatfields 
by minimizing errors in relative photometry of multiply observed stars over the whole survey 
region. Hence, this comparison is an essentially independent test of flatfield corrections 
derived here despite the fact that both catalogs are based on the same observations. We 
bin the photometric differences between the catalogs in 0.01° wide Dec bins and compute 
the median residual for each bin and band. The median value of these medians represent 
zeropoint offsets in each band and are equal to -7, 2, -1, 2 and 5 millimag in the ugriz bands, 
respectively. The rms scatter of the median residuals reflects systematic errors due to flatfield 
errors and we measure 27, 6, 5, 6, and 8 millimag in the ugriz bands, respectively. Except 
in the u band, these values indicate that systematic flatfield errors are very small. In the 
u band, Padmanabhan et al. (2007) expect errors of about 0.01 mag, and the distribution 
width of the s color implies about 0.01 mag for the catalog discussed here, predicting about 
14 millimag instead of the measured 27 millimag. It is plausible that the u band photometry 
may contain systematic errors unrecognized by any of the methods discussed here. 



2.7.2. The tests of catalog quality based on galaxy colors 

The color distribution of galaxies is bimodal (Strateva et al. 2001, Yasuda et al. 2001, 
Baldry et al. 2003). Red galaxies have an especially tight color- redshift relation (Eisenstein 
et al. 2001), with an rms of 0.12 mag for the u — g color, 0.05 mag for g — r, and 0.03 mag 
for the r — i and i — z colors (using model magnitudes). Deviations from the mean relations 
can thus be used to track local calibration problems. Of course, since this is a color-based 
method, it can only constrain d~ mr , and, because red galaxies are faint in the u band, cannot 
achieve high spatial resolution in this band. Nevertheless, it is a useful addition to the 
stellar locus method because it is independent of the Milky Way structure and secondary 
star network (although it is sensitive to errors in the ISM dust extinction correction). 

We select 19,377 red galaxies with SDSS spectra from the redshift range 0.02-0.36 using 
an empirical condition 

< (g - r) - 0.6 - 2.75 x redshift < 0.3, (9) 

and determine their median colors as a function of redshift using 0.01 wide redshift bins. The 
residuals from the median color-redshift relation are then binned by declination to constrain 
SH and by right ascension to constrain <5^f*. The rms for color residuals and the widths 
of distributions of residuals normalized by statistical noise (based on quoted photometric 
errors) are listed in Table 6. 
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The residuals binned in the declination direction are generally small and consistent with 
statistical noise. The largest deviations from are seen for the i — z color, with an rms of 
9 millimag and maximum deviation of 17 millimag (see top panel in Figure [P21 . Although 
the rms is fairly small, it is a factor of 2.9 larger than the expected noise. The shape of 
the i — z residuals is similar to the % — z residuals for stars discussed in § 12.5.21 and shown 
in the bottom panel in Figure [7J The rms scatter for the difference between the stellar and 
galaxy i — z residuals, shown in the bottom panel in Figure [12], is 6 millimag. Since the 
rms for galaxy residuals is larger (9 millimag), it is plausible that the same systematic effect 
dominates the remaining photometric errors for both stars and galaxies. While it is not clear 
what the cause of this (small) effect is, a plausible explanation is the dependence of flatfields 
on source coloi@. 

The results for binning in the right ascension direction are shown in Figure [121 where 
we compare different methods. The rms for implied color errors (with 2° wide bins) from 
galaxies is 0.006 mag for the r — i and i — z colors, 0.012 mag for the g — r color, and 0.018 
mag for the u — g color. The overall behavior of red galaxy color residuals agrees better with 
the PT method, than with color errors implied by the stellar locus method. In particular, 
the large errors outside the —25° < RA < 40° range implied by the latter method are not 
consistent with red galaxy color residuals. On the other hand, both red galaxy colors and 
stellar locus seem to show a trend that the colors are redder around RA=— 10° than around 
RA=40°. The amplitude of this effect varies from about 0.02 mag for the g — r color to 
about 0.01 mag for the i — z color, while the upper limit on such a slope implied by the PT 
comparison is <0.01 mag. 

It is not clear what the cause of this discrepancy is. The obvious culprit is the correction 
for interstellar dust extinction, but the implied deviation is too large to be explained by any 
plausible errors in the SFD98 maps. As shown in Figure [TTJ the median extinction in the r 
band for the — 10 < RA < 40 range is below 0.1 mag, and the resulting median correction for, 
e.g., the g — r color is below 0.04 mag. Hence, to induce a 0.02 mag trend in the g — r color, 
the SFD value for the r band extinction, A r , would have to be in error by 0.05 mag (the 
difference between the values provided for RA=— 10° and RA=40°). This implies relative 
errors for the SFD map in the range from 50% (if A r at RA=— 10° is underestimated) to 
100% (if A r at RA=40° is overestimated), which seems unlikely (although not impossible). 

We conclude that the PT comparison provides a good estimate of the remaining ze- 



19 It is fair to ask whether the applied flatfield corrections, derived from stellar colors, are actually appro- 
priate for galaxies. They are since the i — z color residuals for galaxies without any flatfield corrections are 
about twice as large than those shown in the top panel in Figure [12] 
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ropoint errors in the catalog, as listed in Table 5, but caution that we do not understand 
the above systematic behavior of stellar and galaxy colors, and that only the PT constrains 
possible gray errors. In the next Section, we discuss a comparison to an external dataset, 
that supports this conclusion. 

2.7.3. An external test of catalog quality based on Stetson's standards 

The only large external dataset with sufficient overlap, depth and accuracy to test the 
quality of the Stripe 82 catalog is that provided by Stetson (2000, 2005). Stetson's catalog 
lists photometry in the BVRI bands (Stetson's photometry is tied to Landolt's standards) 
for ~1,200 stars in common (most have V < 19.5). We synthesize the BVRI photometry 
from SDSS gri measurements using photometric transformations of the following form 

mstctson - msoss = Ac i + Bc 2 + Cc + D, (10) 

where mstetson — (BVRI) and m SD ss = {g, 9, r, i), respectively, and the color c is measured by 
SDSS (g—r for the B and V transformations, and r — i for the R and / transformations). The 
measurements are not corrected for the ISM reddening. Traditionally, such transformations 
are assumed to be linear in coloio We use higher-order terms in eq. [TD] because at the 1-2% 
level there are easily detectable deviations from linearity for all color choices (for details 
and plots see Ivezic et al. 2007). The best-fit coefficients for the transformation of SDSS 
gri measurements to the BVRI systerrjEl and low-order statistics for the mstetson — ^sdss 
difference distribution^] are listed in Table 7. We find no trends as a function of magnitude 
at the < 0.005 mag level. With the listed transformations, the SDSS catalog described here 
could also be used to calibrate the data to the BVRI system with a negligible loss of accuracy 
due to transformations between the two systems. 

The BVRI photometry from Stetson and that synthesized from SDSS agree at the level 
of 0.02 mag (rms scatter for the magnitude differences of individual stars; note that the 
systems are tied to each other to within a few millimags by the transformations listed in 
Table 7). This scatter is consistent with the claimed accuracy of both catalogs (the magnitude 



20 For various photometric transformations between the SDSS and other systems, see Abazajian et al. 
(2005) and http://www.sdss.org/dr5/algorithms/sdssUBVRITransform.html, 

21 The same transformations can be readily used to transform measurements in the BVRI system to the 
corresponding gri values because B — V = f(g — r) and R — I = f(r — i) are monotonic functions. 

22 Note that these transformations are valid only for main sequence stars with colors in the range g — r > 0.2 
and r — i < 1.5 (roughly, 0.3< B — V <1.6). Extrapolation outside this range may result in large errors! 
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differences normalized by the implied error bars are well described by Gaussians with widths 
in the range 0.7-0.8). This small scatter allows us to test for the spatial variation of zeropoints 
between the two datasets, despite the relatively small number of stars in common. 

Stars in common are found in four isolated regions that coincide with historical and 
well-known Kapteyn Selected Areas 113, 92, 95, and 113. We determine the zeropoint off- 
sets between the SDSS and Stetson's photometry for each region separately by synthesizing 
BVRI magnitudes from SDSS gri photometry, and comparing them to Stetson's measure- 
ments. The implied zeropoint errors (which, of course, can be due to either the SDSS or 
the Stetson dataset, or both) are listed in Table 8. For regions 1-3 the implied errors are 
only a few millimags (except for the B — g color in region 1). The discrepancies are much 
larger for the three red colors in region 4. A comparison with the results of internal SDSS 
tests described in §2.7.11 and 12.7.21 suggests that these discrepancies are more likely due to 
zeropoint offsets in Stetson's photometry for this particular region, than to problems with 
SDSS photometry. We contacted P. Stetson who confirmed that his observing logs were 
consistent with this conclusion. Only a small fraction of stars from Stetson's list are found 
in this region. 

Given the results presented in this Section, we conclude^! that the rms for the spatial 
variation of zeropoints in the SDSS Stripe 82 catalog is below 0.01 mag in the gri bands. 

3. The Utility of the SDSS Stripe 82 Standard Star Catalog 

As examples of the use of the standard star catalog, we discuss calibration of data 
obtained in non-photometric conditions, and a detailed and robust measurement of the 
morphology of the stellar locus in color-color space. 

3.1. Calibration of Non-photometric Data 

The existence of a technique to photometrically calibrate non-photometric data would 
greatly increase the efficiency of telescopes. As one particular example of how our catalog 
can support a large project, consider the SDSS-II Supernova survey (Sako et al. 2005). This 
survey aims to obtain repeat images of stripe 82 with a sufficient cadence to enable discov- 
ery of new type la supernovae. This requirement sometimes results in observations obtained 



23 Here we assumed that it is a priori unlikely that the SDSS and Stetson's zeropoint errors are spatially 
correlated. 
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through clouds with several magnitudes of extinction. In such highly non-photometric con- 
ditions the standard photometric calibration described in section 12.21 fails because the fields 
with standard stars are too sparsely distributed to be able to resolve fast variations in cloud 
extinction. 



3.1.1. A method to track fast cloud extinction variations 

Due to its high stellar density, the standard star catalog described in this paper can 
be used for calibration of data obtained in grossly non-photometric conditions. The typical 
number of calibration stars in each SDSS field (9x13 arcmin 2 ) at high Galactic latitudes is 
10-15 in the u band, 40-50 in the gri bands, and 30-40 in the z band. Based on tests of 
several non-photometric SDSS-II runs, it was found that the cloud extinction variations can 
be tracked with a sufficient temporal resolution (~3 sec) to obtain photometric zeropoint 
accuracy comparable to that characteristic for photometric nights (1-2% in gri and 2-3% in 
u and z; Ivezic et al. 2004a). 

The calibration is done in two steps. First, the implied zeropoints (whose variation is 
dominated by cloud extinction), zp, defined by 

m ca ii bra ted = -2.5 * \og(counts) + zp, (11) 

are computed for each star and median filtered in time using a window with 5 stars in 
order to avoid outliers. Note that we assume that clouds are gray and do not allow for 
color terms, an assumption which is justified a posteriori (see §3.1.4p . In the second step, 
zeropoints are evaluated for each 2048 pixel wide (cross-scan direction) and 100 pixel long 
(in-scan direction) image segment, hereafter called a calibration patch (not to be confused 
with secondary star patches discussed in §2.2p . That is, a calibration patch is a ~9 arcmin 2 
large rectangle with an aspect ratio 1:20, and the zeropoints are evaluated every 2.6 seconds 
of time (but note that the variations are smoothed out by the 54 sec long exposure time). 

The patch is much narrower in the in-scan direction because tests have shown that 
zeropoint gradients across a field are much larger, by a factor 10-50, in this direction (see 
Figure [T4"]) . Consider three stars, star A, a star B that is, say, 25 arcmin (the column-to- 
column separation) away from the star A in the scan direction, and a star C that is 25 arcmin 
away from the star A in the cross-scan direction. Stars A and C are observed at the same 
time and the difference in their implied zeropoints measures the structure function of cloud 
opacity on a 25 arcmin spatial scale. This is true irrespective of the cloud motion relative to 
the boresight. Here, the structure function of cloud opacity (SF hereafter) is defined as the 
rms width of the distribution of zeropoint differences evaluated for pairs of points separated 
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by some distance. 

On the other hand, stars A and B are observed at times that differ by 1.7 minutes. If 
the component of the cloud angular velocity on the sky relative to the boresight and parallel 
to the scanning direction is u ° min -1 , the zeropoint difference for stars A and B samples the 
cloud structure on spatial scales of 25 uj/uo s arcmin, where u; s =0.25 ° min -1 is the sidereal 
scanning rate (here for simplicity we assumed uj 3> u s , which is supported by the data). The 
observed behavior of zeropoints, such as that shown in Figure [H], implies wind velocity in 
the rang uj=3 — 15 min 1 , or uj/uj s ~ 12 — 60. Hence, drift scanning has the unfortunate 
property that the motion of an inhomogeneous extinction screen with a speed much larger 
than the sidereal scanning rate greatly magnifies the effective zeropoint variations in the 
scan direction. 

The zeropoints for each calibration patch are computed by taking all the stars from 
the patch, or finding the closest three stars for sparsely populated patches, and adopting 
the median value of their zeropoints. This is certainly not the only, nor perhaps the ideal 
approach to calibrate patches, but we found that it works well in practice. The zeropoint 
error is evaluated from the root-mean-square scatter of zp evaluated for each calibration star, 
divided by the square root of the number of stars. We now discuss the performance of this 
method. 



3.1.2. Performance and Quality Tests 

The top panel in Figure [T5] summarizes the behavior of cloud extinction in the r band, 
as measured by the zeropoint zp discussed above, for an SDSS-II SN run (5646) obtained 
in strongly non-photometric conditions. Although the cloud extinction during the first 90 
minutes (corresponding to 150 SDSS fields) varies between and ~6 mag, it is possible to 
robustly calibrate these data. Figure [TBI zooms in on a 8 minute stretch of the same data 
where the cloud extinction varies between and ~3 mag, with changes as fast as 0.05 mag 
s -1 (almost 2 mag per SDSS field!). As shown in the figure (middle left panel), the residuals 
have a distribution width of only 0.07 mag. The middle right panels in Figures [15] and 
[TBI demonstrate that most of this scatter is contributed by random photometric errors (i.e. 
errors in extracted source counts), rather than by calibration errors (large cloud extinction 
results in a smaller number of calibration stars, as well as in a lower SNR for those calibration 



24 This range is equivalent to angular speeds of up to a half of the Moon's diameter per second. The 
plausibility of this wind velocity range was verified in extensive visual observations of the full Moon during 
frequent grossly non-photometric nights in Seattle. 
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stars that are detected). Even with such a large and rapidly varying cloud extinction, the 
zeropoint errors are smaller than 0.05 mag, with a median value of less than 0.02 mag. An 
example of a run with somewhat thinner and much more stable cloud cover is shown in 
Figure [T71 

The calibration performance in other bands is similar. For example, although the num- 
ber of calibration stars is smaller in the u band than in the r band, the median zeropoint 
error for the same stretch of data as shown in Figure [T71 is still only 0.01 mag, as illustrated 
in Figure [T51 

3.1.3. The Summary of Calibration Accuracy 

A summary of the final zeropoint errors as a function of cloud extinction and band 
for one of the worst runs is shown in the left column in Figure [T^J As the figure shows, 
the data can be calibrated with small zeropoint errors even for such a bad case. Typically, 
the zeropoint errors, for the same cloud extinction, are about twice as small as in this run. 
A calibration summary for a run with optically thick, but exceptionally smooth, clouds is 
shown in the right column in Figure [191 Overall, for cloud extinction of X mag, the zeropoint 
uncertainty is typically smaller than (0.02 — 0.05)X for 95% of calibration patches, with a 
median of (0.01-0.02)X . 

3.1.4- The Cloud Color and Structure Function 

We detect no dependence of the calibration residuals on the stellar color or cloud thick- 
ness at a few millimag level. This is consistent with the lack of selective extinction by clouds. 
The lack of a color correlation in the u and z bands implies that the well-known cloud grey- 
ness extends beyond visual wavelengths. Another method to quantify selective extinction by 
clouds is to directly compare zeropoints from different bands. As shown in Figure [201 the 
cloud extinction is similar in all bands for most fields. A few cases where there are deviations 
of a few tenths of a magnitude can be easily understood as due to temporal changes in the 
cloud opacity (recall that the data from different bands are obtained over ~5 minutes of 
time). 

The calibration accuracy is determined by the size of the calibration patches. For ex- 
ample, a smaller patch would suffer less from the spatial variation of cloud extinction, but 
it wouldn't have enough stars to beat down the noise of their individual photometric mea- 
surements (~0.02 mag for sufficiently bright stars). The detailed scaling of this accuracy 
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with patch size depends on the cloud spatial structure function. The geometry of the SDSS 
camera allows us to study the cloud SF on scales exceeding 2°. Figure [21] compares zero- 
points in different columns for two runs with significantly different cloud behavior. While 
zeropoints from different columns generally track each other, there can be differences exceed- 
ing a magnitude (they generally scale with the cloud optical thickness). These differences 
increase with the distance between the camera columns. Figure [221 shows a typical behavior: 
for small spatial scales (<2°) the SF is roughly a linear function of distance, and it scales 
roughly linearly with the cloud extinction. At a 1° scale, the SF is typically of order 2-10% 
of the cloud extinction. For example, even for clouds 3 mag thick, the SF at 2 arcmin scales 
is typically <0.01 mag. 



3.1.5. Implications for Surveys such as LSST 

The Large Synoptic Survey Telescope (LSST) is a proposed imaging survey that will 
attempt to maximize its observing time by accepting non-photometric conditions. At the 
same time, it has adopted exquisite requirements for its photometric accuracy, including 1% 
accuracy of its internal photometric zeropoint errors across the sky. Our analysis allows 
us to answer the following question: "What is the largest cloud extinction that still allows 
photometric calibration accurate to 1%?". 

A similar approach to the calibration of LSST data as presented here (assuming that a 
standard star catalog is available, e.g. from prior demonstrably photometric nights) would 
benefit from several effects: 

1. The LSST will not use a drift-scanning technique and thus the calibration patches can 
be squares; for the same area, this results in a ~5 times smaller angular scale (~ 3 
arcmin), compared to the 1:20 rectangles we have used to calibrate SDSS drift-scanning 
data. On these angular scales, the cloud structure function is roughly linear and thus 
the zeropoint error would be ~5 times smaller, or of the order 1% or less through 
clouds as thick as 1 mag (conservatively assuming that SDSS errors would be 0.05A, 
see § I3.1.3p . We note that the shorter exposure time for LSST (30 seconds, or about 
a factor of two shorter than for SDSS) is not relevant because clouds would typically 
move by more than a degree during the exposure. This is more than an order of 
magnitude larger distance than the size of the calibration patch and thus the structure 
function analysis remains valid. 

2. LSST data will be deeper than SDSS by about 2-3 mag. With a conservative assump- 
tion that logiV oc 0.3m for faint stars (0.6 for Euclidean counts), the surface density 
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of calibration stars will be about ten times larger for LSST than for SDSS. This larger 
density enables ten times smaller patches, or about three times smaller angular scale for 
calibration (~ 1 arcmin), resulting in another factor of three improvement of accuracy. 

3. Fitting a smooth function for cloud opacity over several calibration patches would 
result in further improvements. 



The first two points predict that LSST data could be calibrated with the required 1% ac- 
curacy even through 3 mag thick clouds. Given the various extrapolations, we conservatively 
suggest the range of 1-3 mag as the upper limit on the acceptable cloud opacitjQ 



3.2. The Morphology of the Stellar Locus 

The improved accuracy of averaged photometry provides " crisper" color-color diagrams 
and also reveals new morphological features. An example of such a color-color diagram is 
shown in Figure |2"3"1 

This is a similar plot to Figure 1 from Smolcic et al. (2004), except that only non- 
variable point sources are shown (note the absence of quasars) and averaged photometry is 
used. The white dwarf/M dwarf "bridge" discussed by Smolcic et al. is clearly visible, as 
well as the locus of probable solar metallicity giants ( this identification is based on models, 
e.g. Kurucz 1979) which bifurcates from the main locus at u — g ~ 2.5 and g — r ~ 1. Note 
also the well-defined BHB locus (u — g ~ 1.1 and g — r from —0.3 to 0.1) and the white 
dwarf locus (u — g ~ 0.35 and g — r from —0.3 to ~0.0). A new locus-like feature, that is 
not visible in Figure 1 from Smolcic et al., is discernible at u — g ~ and g — r ~ —0.2. The 
great value of the accurate u band photometry is clearly evident at e.g. g — r = —0.2: the 
u — g color distribution is tri-modal! The bluest branch is consistent with He white dwarfs, 
and the middle branch with hydrogen white dwarfs, as supported by Bergeron et al. (1995) 
white dwarf models and detailed analysis of SDSS spectra (Eisenstein et al. 2006). The 
reddest branch is made of blue horizontal branch stars (see Sirko et al. 2004 and references 
therein) . 

The exciting fact that one can distinguish He and H white dwarfs using photometry 
alone is a consequence of the improved photometric accuracy due to averaging many epochs. 
Figure [24] reiterates that point. Note the striking difference between the two bottom panels: 



25 0f course, cloud opacity decreases the imaging depth and data with clouds thicker than ~1 mag may 
be undesireable for reasons other than calibration accuracy. 
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while one could be convinced that the He white dwarf sequence is a real feature in the 
bottom right panel, its existence is clearly evident when using the improved photometry, as 
shown in the bottom left panel. In summary, the multi- epoch observations provide both the 
identification of variable sources and much more accurate colors for non-variable sources. 
This bodes well for science deliverables from upcoming large-scale imaging surveys. For 
example, LSST will obtain over its 10-year long mission similar repeat imaging as discussed 
here, but about 2.5 mag deeper, with about 100 or more observations per band and object, 
and over about two orders of magnitude larger area. Although these new surveys will not 
have a spectroscopic component like SDSS did, the multi-epoch nature of their imaging will 
provide alternative information-rich datasets. 



4. Discussion and Conclusions 

Using repeated SDSS measurements, we have constructed a catalog of over a million 
candidate standard stars. The catalog is publicly available from the SDSS websita^l. Several 
independent tests suggest that both internal zeropoint errors and random photometric errors 
for stars brighter than (19.5, 20.5, 20.5, 20, 18.5) in ugriz, respectively, are at or below 
0.01 mag (about 2-3 times better than for an individual SDSS scan). This is by far the 
largest existing catalog with multi-band optical photometry accurate to ~1%, and breaks 
the accuracy barrier discussed by e.g. Stubbs & Tonry (2006, and references therein). These 
observations were not obtained for the specific purpose of calibration, but were part of the 
regular SDSS observational program. When compared to, for example, the heroic calibration 
efforts by Landolt, Stetson, and others, it seems justified to call the method presented here 
"industrial" photometry. However, the catalog presented here is not without its problems. 

The selection of candidate stars was simply based on the absence of variability. It is fairly 
certain that not all variable sources are recognized because of the limited number of repeated 
observations (~10). For example, an eclipsing binary with much shorter eclipse duration than 
the orbital period could easily escape detection. Furthermore, some of these sources may 
not even be stars. A cross-correlation with the SDSS spectroscopic database yields 99,000 
matches in the overlapping region, including 44,000 spectra classified as stars. About 70 
candidate standard stars are actually spectroscopically confirmed quasars! Apparently, a 
small fraction of quasars (a few percent, for a detailed analysis see Sesar et al. 2007, in 
prep.) cannot be detected by variability (at least not using the number of epochs, their 
time distribution, and photometric accuracy employed in this work). Indeed, we have also 



26 See 



http: / / www.sdss.org/dr5 / products / value ..added/index. html 
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found three spectroscopically confirmed SDSS quasars among Stetson's standards which 
were observed 20-30 times and which showed no variation. Similarly, about 300 candidate 
standard stars have SDSS spectra classified as galaxies. Nevertheless, the inspection of color- 
color diagrams strongly suggest that the overwhelming majority of the standard stars are 
found on the stellar locus. 

Remaining systematic errors are another important concern. Effectively, we have as- 
sumed that PT problems average out in many patches when deriving flatfield corrections 
using stellar colors. This may not be true at a level not much smaller than 1%, and thus 
the remaining gray problems at such a level may be present in the catalog. Despite these 
residual problems, we believe that internal consistency of the catalog (i.e. when ignoring A m 
from eq. H|) is such that the rms width for the function 5 m (RA, Dec) from eq. H] evaluated for 
all stars in the catalog is at most 0.01 mag in the griz bands and perhaps just slightly larger 
in the u band (very unlikely exceeding 0.02 mag). In addition to gray problems and overall 
flatfield errors, the dependence of flatfields on source color is probably the largest remaining 
systematic error. 

We illustrate several use cases for this catalog, including the calibration of highly non- 
photometric data and robust selection of stars with peculiar colors. We find that LSST 
and similar surveys will be able to observe in partially cloudy (non-photometric) nights 
because even cloudy data can be accurately calibrated with a sufficiently dense network of 
calibration stars. Such a dense network will be self-calibrated by LSST very soon after first 
light, using an approach developed for SDSS data by Padmanabhan et al. (2007). Given 
such a network, SDSS experience suggests that LSST can maintain its required photometric 
calibration accuracy of 1% even when observing through 1-3 mag thick clouds. 

Perhaps the most exciting conclusion of this work is that it may become obsolete in only 
a few years due to the advent of next-generation surveys such as Pan-STARRS and LSST. 
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Table 1. The Color-term Corrections" 



column 


(du) 


0~du 


(dg) 


a dg 


(dr) 




(di) 




(dz) 




1 


1.0 


0.3 


1.8 


0.4 


-1.2 


0.3 


0.0 


0.2 


0.6 


0.5 


2 


11.0 


2.3 


7.6 


1.6 


-1.1 


0.2 


-0.0 


0.2 


-0.8 


1.3 


3 


-5.4 


0.8 


-2.0 


0.0 


2.9 


1.0 


-0.2 


0.3 


-0.2 


0.5 


4 


-9.5 


2.1 


-1.0 


0.2 


-0.9 


0.1 


0.0 


0.0 


-0.4 


0.3 


5 


3.9 


1.0 


-4.4 


1.3 


-0.0 


0.2 


-0.2 


0.1 


0.4 


0.7 


6 


-0.5 


0.8 


-2.6 


0.4 


-0.6 


0.1 


-0.2 


0.3 


0.9 


1.2 



a The median and rms scatter for photometric corrections that place the 
measurements in six camera columns on the survey system (in millimag). 



Table 2. The s,w,x,y Principal Color Definitions 



PC 


A 


B 


C 


D 


E 


F 


s 


-0.249 


0.794 


-0.555 


0.0 


0.0 


0.234 


w 


0.0 


-0.227 


0.792 


-0.567 


0.0 


0.050 


X 


0.0 


0.707 


-0.707 


0.0 


0.0 


-0.988 


y 


0.0 


0.0 


-0.270 


0.800 


-0.534 


0.054 
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Table 3. The Effect of Repeated Measurements on the Width of the Stellar 

Locus 



PC 


rms 1 obs a 


rms N obs 6 


median c 


width for PC/PCerror d 


s 


31 


19 


3.0 


3.02 


w 


25 


10 


1.1 


1.82 


X 


42 


34 


1.2 


5.34 


y 


23 


9 


0.8 


1.64 



a The locus width determined using single-epoch SDSS observations (in mil- 
limag) . 

b The locus width determined using multiple SDSS observations (in mil- 
limag) . 

c The median principal color determined using multiple SDSS observations 
(in millimag). 

d The locus width normalized by expected measurement errors. 



Table 4. The Flatfield Corrections" 



band 


width 6 


min c 


max d 


ur 


22 


-53 


53 


gr 


12 


-27 


19 


r 


7 


-17 


17 


ri 


4 


-10 


13 


rz 


7 


-14 


19 



a The r band correction is determined using observations by the Photomet- 
ric Telescope, and the ugiz corrections are determined using the stellar locus 
method (see 32XT1 and 12X21) . 

b The root-mean-square scatter for applied flatfield corrections (in millimag). 

c The minimum value of the applied correction (in millimag). 

d The maximum value of the applied correction (in millimag). 
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Table 5. The statistics of the median PT-2.5m residuals 



band 


<> a 


width 6 


min c 


max d 


N e 


<>f 


u 


-2 


9 


-18 


23 


175 


5 


9 


6 


7 


-4 


17 


647 


5 


r 


3 


7 


-7 


10 


627 


3 


i 


4 


7 


-10 


17 


621 


2 


z 


1 


8 


-16 


15 


286 


-2 



a The median value for the bin medians (in millimag). There are 24 bins, 
distributed inhomogeneously in the right ascension direction. 

b The root-mean-square scatter for the bin medians (in millimag). 

c The minimum value for the median residuals (in millimag). 

d The maximum value for the median residuals (in millimag). 

c The median number of stars per bin. 

f The median value of the residuals for stars with colors within 0.02 mag from 
the crossing colors (in millimag). 
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Table 6. Residuals from the mean color- redshift relation 61 for red galaxies. 



color 


rms(Dec) b 


x(Dec) 


rms(RA) c 


x(RA) 


u-g 


21 


1.5 


18 


1.3 


g-r 


4 


1.3 


12 


2.6 


r — i 


3 


1.3 


6 


3.5 


i — z 


9 


2.9 


6 


2.9 



a The table lists the rms widths of color-residual distributions (in millimag), 
and the widths of distributions of residuals normalized by statistical noise (x), 
using mean color- redshift relations (see text). 

b The rms for the declination direction, using 0.025° wide bins (0.1° for the 
u — g color). 

c The rms for the right ascension direction, using 2° wide bins (5° for the 
u — g color). 
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Table 7. Comparison with Stetson's standards: I. griz to BVRI transformations 



color 


<-> b 


u med 




<> all 


a all 


k9 


B9 






B-g 


-1.6 


8.7 


1.4 


1.0 


32 


0.2628 


-0.7952 


1.0544 


0.0268 


V-g 


0.8 


3.9 


1.0 


0.9 


18 


0.0688 


-0.2056 


-0.3838 


-0.0534 


R-r 


-0.1 


5.8 


0.9 


1.2 


15 


-0.0107 


0.0050 


-0.2689 


-0.1540 


I -i 


0.9 


6.1 


1.0 


1.2 


19 


-0.0307 


0.1163 


-0.3341 


-0.3584 



a These transformations are valid only for main sequence stars with colors in the range g — r > 0.2 
and r — i < 1.5 (roughly, 0.3< B — V <1.6). Extrapolation outside this range may result in large 
(0.1 mag or larger for hot white dwarfs) errors! 

b The median value of median transformation residuals (differences between the measured values 
of colors listed in the first column and those synthesized using eq. [TOl in 0.1 mag wide g — r bins 
for stars with 0.25< g — r <1.45 (in millimag). These medians of medians measure the typical 
level of systematics in the gri-to-BVRI photometric transformations introduced by the adopted 
analytic form (see eq. [TO]) . 

c The root- mean-square scatter for median residuals described above (in millimag). 

d The root-mean-square scatter for residuals normalized by statistical noise. The listed values 
are ~1, which indicates that the scatter around adopted photometric transformations listed under 
b) is consistent with expected noise. 

c The median value of residuals evaluated for all stars (in millimag). 

f The root-mean-square scatter for residuals evaluated for all stars (in millimag). 

Coefficients A-D needed to transform SDSS photometry to the BVRI system (see eq. [TOl . 
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Table 8. Comparison with Stetson's standards: II. Photometric zeropoint variations 



color 


<>m 






<>m 


<2 




<>R3 


<3 


N C R3 


<>m 






B-g 


-29 


21 


92 


6 


27 


165 


8 


42 


155 


-4 


27 


281 


V-g 





17 


99 





15 


217 


6 


25 


161 


17 


19 


282 


R-r 


-6 


16 


58 


4 


16 


135 


-8 


12 


11 


39 


27 


60 


I -i 


-11 


16 


94 


6 


18 


205 


2 


16 


124 


19 


15 


47 



a The median value of residuals (in millimagnitudes) for transformations listed in the first column, 
evaluated separately for regions 1-4, defined as: Rl: RA~325, Dec<0; R2: RA~15; R3: RA~55; 
R4: RA~325, Dec>0. 

b The root-mean-square scatter for the transformation residuals (in millimagnitudes). 

c The number of stars in each region with good photometry in the required bands. 
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Fig. 1. — The median magnitude error as a function of magnitude (left column, symbols), the x 2 
per degree of freedom distribution (solid line, middle column), and the number of observations in 
each band (right column) for candidate standard stars from SDSS stripe 82 (selected by x 2 < 3 
in the gri bands). The solid lines in the left column show cumulative magnitude distribution 
normalized to 10 at the faint end. The dashed lines in the middle column show the x 2 P er degree 
of freedom distribution for a Gaussian error distribution and nine degrees of freedom. Its similarity 
with the measured distributions suggests that the magnitude errors computed by the photometric 
pipeline are reliable (they may be slightly underestimated, by about 10%, in the u band). 



-41 - 




Fig. 2. — A comparison of the color-magnitude and color-color distributions for variable and 
non-variable unresolved sources from SDSS stripe 82. The distributions of sources with a root- 
mean-square (rms) scatter of the g band magnitude below 0.05 mag are shown by linearly-spaced 
contours (non- variable sources). Sources with 320° < RA < 330°, with rms scatter in the g band 
larger than 0.05 mag, and x 2 greater than three, are shown by dots (variable sources). The dots 
are color-coded according to the observed rms scatter in the g band (0.05-0.10 mag, see the legend, 
red indicates variability larger than 0.1 mag). Note how low-redshift quasars (u — g < 0.6) and RR 
Lyrae (u — g ~ 1.1, g — r ~ 0; Ivezic et al. 2005) clearly stand out as variable sources (red points). 
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Fig. 3. — The color-magnitude and color-color distributions of candidate standard stars 
brighter than g = 20. The dots are color-coded according to the observed rms scatter in the 
g band (0-0.05 mag, see the legend). Note the absence of quasars and RR Lyrae visible in 
Figure [21 
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Fig. 4. — Analogous to Figure [31 except that only variable sources are shown, and with a 
different color coding (0.05-0.10 mag range of the g band rms, instead of 0-0.05 mag). 
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Fig. 5. — The top panel shows the sky density of all the point sources with at least 4 
observations in each of the g, r and i bands (solid), and of selected candidate standard (non- 
variable) stars (dashed). The bottom panel shows the ratio of the two curves shown in the 
top panel. For reference, Galactic coordinates, (l,b), are (46,-24), (96,-60) and (190,-37) 
for a J2000 =-50°, 0° and 60° (at 5 J20 oo=0 ). 
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Fig. 6. — The top panel compares the five transmission curves, (f) (see eq]3]for definition), that 
define the SDSS photometric system (dashed lines), to measured transmission curves sorted 
by camera column and shown by solid lines, as marked in the panel (the area under each 
curve is unity by definition). The bottom panel shows the differences between the measured 
transmission curves and the curves that define the system. These differences induce color 
terms that result in systematic photometric errors as a function of source color. The largest 
color terms are present for the z band in camera columns 2 and 6, with the errors well 
described by Az 2 = z 2 — ^sdss = — 0.019(r — i) and Az 6 = +0.017(r — i) mag, respectively. 
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Fig. 7. — A test based on the stellar locus position angle in the i — z vs. r — i color-color 
diagrams that demonstrates the existence of color terms between different camera columns. 
The y axis shows the difference in the i — z color residuals of blue (0.1< r — i < 0.2) and 
red (0.8< r — i <1.4) stars, with the residuals computed as the difference between measured 
i — z colors and those predicted using the mean stellar locus (see the bottom right panel 
in Figure [3]). The vertical solid lines mark the approximate boundaries between different 
camera columns, with the vertical dashed lines marking approximate boundaries between 
the "north" and "south" strips in a stripe (see §2.11 for definitions). The top panel shows 
results before applying corrections for different transmission curves (see Figure EJ) and the 
bottom panel shows results based on corrected photometry. As evident, the residuals in the 
bottom panel are much smaller, with rms scatter decreasing from 9 millimag to 3 millimag. 
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Fig. 8. — The top two panels show two projections of the stellar locus (rotated locus from the 
ugr and gri planes, see Figure [3]) constructed using averaged photometry. The large green dots 
show the median values of the s and w principal colors (perpendicular to the locus at its blue edge) 
in bins of the principal color along the locus. The middle and bottom rows show histograms for 
each principal color on linear and logarithmic scales (essentially the locus cross-sections). The blue 
(narrow) histograms are constructed using the averaged photometry, and the magenta histograms 
are based on single-epoch photometry. The former are narrower than the latter due to increased 
photometric accuracy. The best-fit Gaussians, with parameters listed in the middle row, are shown 
by the dashed lines in the bottom row (only for the averaged photometry). 
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Fig. 9. — The top panel shows the distribution of the residuals between the PT photom- 
etry and averaged magnitudes in the gri bands. The second panel shows applied flatfield 
correction in the r band, which was derived as the mean of the residuals shown in the top 
panel. The remaining two panels show the applied flatfield corrections in the other four 
bands, expressed relative to the r band (middle: u, bottom: gri), which were derived using 
the stellar locus colors. The low-order statistics for these corrections are listed in Table 4. 
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Fig. 10. — The panels show the dependence of the position of the stellar locus in ugriz color 
space, as parametrized by the median principal colors swxy, as a function of RA. Close to the 
edges, the median colors deviate significantly from 0. This is caused by intrinsic changes in 
the stellar locus due to stellar population variations and overestimated interstellar extinction 
corrections for red stars, rather than calibration problems. 
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Photometric zeropoint errors implied by the stellar locus position 
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Fig. 11. — The top panel shows the implied photometric zeropoint errors based on the stellar locus 
method {ugiz from the bottom to the top at either edge). While the implied errors are small for 
—25 < RA < 40, they become exceedingly large outside this range. This is due to problems with 
the stellar locus method rather than due to problems with calibration. The second panel shows the 
median r band extinction derived from the Schlegel, Finkbeiner & Davis (1998) maps. The bottom 
two panels show the median residuals between the PT photometry and averaged magnitudes in the 
uz (third row) and gri bands (bottom row). The low-order statistics for these residuals are listed 
in Table 5. 
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Fig. 12. — The top panel shows the dependence of the median i — z color residuals (with 
respect to a mean color-redshift relation) for red galaxies as a function of Dec. The vertical 
solid lines mark the approximate boundaries between different camera columns, with the 
vertical dashed lines marking approximate boundaries between the "north" and "south" 
strips in a stripe. The small rms scatter of only 9 millimag demonstrates that flatfield 
corrections based on the stellar locus position in color space are also applicable for galaxies. 
The bottom panel shows the difference between the values shown in the top panel and the 
curve shown in the bottom panel in Figure The rms scatter for the residuals in the bottom 
panel is 6 millimag. 
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u-g errors: stellar locus (line), galaxies (triangles) and PT (circles) 
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Fig. 13. — The comparison of systematic color errors implied by different methods. Note 
that the errors implied by the stellar locus method (line) become very large outside the 
—25° < aj2ooo < 40° range. As galaxy colors (triangles) and a direct comparison with SDSS 
secondary standard star network (circles) suggest, this is due to problems with the stellar 
locus method rather than due to problems with calibration. 
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Fig. 14. — The comparison of cloud extinction gradients in the in-scan (RA, horizontal axis) 
and cross-scan (Dec, vertical axis) directions for SDSS run 5759, on a spatial scale of ~0.4° 
(each point is derived using zeropoints from three calibration patches). Note the different 
axis scales. For this particular run, the distribution width is 9.6 times larger for the in-scan 
than for the cross-scan direction. This is a consequence of cloud motion relative to the 
boresight, and the drift-scanning technique. 
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Fig. 15. — The top panel summarizes the behavior of cloud extinction in the r band over 1.5 hours during 
SDSS-II run 5646. Individual calibration stars are shown by dots, and the adopted zeropoint is shown by the 
line. The calibration residuals for each star are shown by dots in the second panel. The root-mean-square 
scatter for these residuals evaluated for each field is shown by the line. The distribution of the residuals is 
shown in the left panel in the third row as the solid line. The median and equivalent Gaussian a evaluated 
from the inter-quartile range are also shown in the panel, as well as a Gaussian corresponding to these 
parameters (dashed line). The right panel is analogous, except that the residuals are normalized by the 
expected errors. The distribution of implied cloud extinction is shown in the bottom left panel, and the 
distribution of standard errors for the adopted photometric zeropoints (computed from the rms width of the 
distribution of residuals) is shown in the bottom right panel (a few points, about 4 out of 1800 calibration 
patches, are off scale). 
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Fig. 16. — Analogous to Fig. [Jj)J except that only ~8 minutes of data with large cloud 
extinction is shown. Note that the changes in cloud extinction are resolved down to time 
scales well below 1 minute. 
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Fig. 17. — Analogous to Fig. [T51 except that 1.5 hours of data from run 5759, which had 
somewhat thinner and much more stable cloud cover, is shown. Note that the median 
photometric zeropoint error is below 0.01 mag, although the median cloud extinction is 
larger than 1 mag. 
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Fig. 18. — Analogous to Fig. [T7J except that the u band calibration summary is shown. 
Despite the smaller number of calibration stars than in redder bands, and over a magnitude 
of cloud extinction, it is still possible to photometrically calibrate these data with a median 
error of only ~0.01 mag. 
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Fig. 19. — A summary of calibration accuracy as a function of the cloud extinction and 
band (urz as marked in the panels). Each small symbol represents one calibration patch (a 
~9 arcmin 2 large rectangle with 1:20 aspect ratio). Zeropoint error is determined from the 
root-mean-scatter of photometric residuals. The large symbols show the median zeropoint 
error in 1 mag wide bins of cloud extinction. The left column shows data for one of the 
photometrically worst SDSS-II SNe runs (5646), and the right column is for a run with 
optically thick, but exceptionally smooth, clouds (5759). Note that the data can be calibrated 
with zeropoint errors typically smaller than a few percent even through clouds several mag 
thick. 
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Run 5646, column 1, implied cloud color 
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Fig. 20. — The color of cloud extinction in SDSS bands. Each symbol represents one field and shows 
the difference in cloud extinction between the two bands as a function of the r band extinction. 
The measurements in different bands are obtained over ~5 minutes of time and thus even gray 
clouds with spatially varying extinction could produce the observed non-gray (non-zero) values. 
The dashed lines in the bottom two panels indicate the expected correlation if the color variations 
are due to temporal changes in the gray cloud thickness (rather than due to intrinsic color changes). 
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Fig. 21. — A comparison of cloud extinction independently measured for six camera columns. The 
top two panels show the 30 minutes of measurements in the g and r band for the same run (5759) 
shown in Figure [TTJ Individual camera columns are color-coded according to the legend shown in 
the top panel in Figure The third panel shows the difference of the r band zeropoints measured 
in one of the edge columns (6) and zeropoints from other five columns. The bottom two panels are 
analogous to the second and third panel, except that the data are from a run with exceptionally 
patchy clouds (5646, the first 30 minutes of the data from Figure [T5l are shown). Note the varying 
scale for the y axis. 
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Run 5759, cloud SF, r band 
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Fig. 22. — The cloud structure function in the r band for run 5759 and for the same stretch of data 
as shown in Figure [T7l (the median cloud extinction is 1.3 mag). The circles show the rms width 
of the distribution of zeropoint differences between camera column 1 and (five) other columns. 
This width is corrected for a 0.015 mag contribution from the measurement errors and shown by 
squares. The triangles show the width of the distribution of zeropoint differences in the in-scan 
direction, with the distance scale multiplied by 30. This multiplication factor measures the cloud 
speed relative to the boresight in the in-scan direction (see §3.1.1j) . 
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Fig. 23. — The g — r vs. u — g color-color diagrams for all non-variable point sources 
constructed with the improved averaged photometry (dots). Various stellar models (Kurucz 
1979; Bergeron et al. 2005; Smolcic et al. 2006) are shown by lines, as indicated in the 
figure. 
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Fig. 24. — An illustration of the advantages of repeated photometric measurements. The 
top left panel shows the blue corner of the g — r vs. u — g diagram from Figure [23] for all 
point sources with the averaged photometry. The top right panel shows only the variable 
sources (dominated by low-redshift quasars), and the bottom left panel shows the non- 
variable sources (dominated by white dwarfs), classified using low-order lightcurve moments. 
The bottom right panel shows the same non-variable sources, but using their DR5 single- 
epoch photometry. A comparison of the bottom two panels shows the striking improvement 
made possible by the use of multiple observations of the same field. 



